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1 Introduction 

In this paper, we consider the following compressible and isentropic Navier-Stokes equations 
with density-dependent viscosities 

[pt + (pu)x = 0, x e R, t > 0, 

\{pu) t + (pu 2 +p(p)) x = {p{p)u x ) x , 

where p(t, x) > 0, u(t, x) represent the density and the velocity of the gas, respectively. Let the 
pressure and viscosity function be given by 

p(p)=Ap', p(p) = Bp a , (1.2) 

respectively, where 7 > 1 denotes the adiabatic exponent, a > and A, B > are the gas 
constants. For simplicity, it is assumed that A = B = 1. 

Consider the Cauchy problem (jl.ip with the initial values 

(p,pu)(Q,x) = (p ,m )(x) ->■ (p±,m±), as x ->■ ±00, (1.3) 

where p-i- > 0, m-t are prescribed constants. 

The large time behavior of solutions to (jl.ip - (|1.3p is expected to be closely related to that 
of the Riemann problem of the corresponding Euler system 



Pt + (pu) x = 0, 
(pu) t + (pu 2 +p(p)) 3 



;i-4) 



with Riemann initial data 



„ (p_, m_), x < 0, 
( P ,p U )(0,x):=(p5,m5)(x)=<; (1.5) 

(p + ,m+), x > 0. 

Different initial states (jl.5p produce different type of waves, namely, shock waves and rar- 
efaction waves for the one-dimensional compressible isentropic Euler equations (|1.4j) . However, 
as pointed out by Liu-Smoller [25], among the two nonlinear waves, i.e., shocks and rarefaction 
waves, only rarefaction waves can be connected to vacuum. When vacuum appears, the stability 
of rarefaction waves to the ID compressible Navier-Stokes equations is an important issue. 

When the viscosity p(p) is a constant, there have been extensive studies on the stability of the 
rarefaction waves to the ID compressible Navier-Stokes equations under the assumptions that 
the rarefaction waves and the solutions are away from the vacuum (see [8], [T7], [21], [23], [26] . 
|29| . [30] and the references therein). However, when vacuum appears, the well-known results by 
Hoff-Serre [13], Xin [37] and Rozanova |34] show that the solutions of the compressible Navier- 
Stokes equations with constant viscosity may behave singularly, in particular, in the case that 
the fluids jump to far field vacuum. Liu, Xin and Yang first proposed in [27] some models of 
the compressible Navier-Stokes equations with density-dependent viscosities to investigate the 
dynamics of the vacuum. On the other hand, when deriving by Chapman-Enskog expansions 
from the Boltzmann equation, the viscosity of the compressible Navier-Stokes equations depends 
on the temperature and thus on the density for isentropic flows. Also, the viscous Saint- Venant 
system for the shallow water, derived from the incompressible Navier-Stokes equations with 
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a moving free surface, is expressed exactly as in (|l.ip - (jl.2p with a = 1 and 7 = 2 (see [9]). 
However, there appear new mathematical challenges in dealing with such systems. In particular, 
these systems become highly degenerate. The velocity cannot even be defined in the presence 
of vaccum and hence it is difficult to get uniform estimates for the velocity near vacuum. The 
global existence of generak weak solutions to the compressible Navier-Stokes equations with 
density-dependent viscosities or the viscous Saint- Venant system for the shallow water model 
in the multi-dimensional case remains open, and one can refere to [4], [5], |10j . |31j for recent 
developments along this line. 

There are a large number of literatures on mathematical studies of f)l.l[) - f)1.2|) with various 
initial and boundary conditions. If the initial density is assumed to be connected to vacuum 
with discontinuities, Liu, Xin and Yang first obtained in [27] the local well-posedness of weak 
solutions. The global well-posedness was obtained later by [E], p2], [33], [38] respectively. 
The case of initial densities connecting to vacuum continuously was studied by [7J, [36] . |38| 
and [39] respectively. However, most of these results concern with free boundary problems. 
Recently, initial-boundary- value problems for the one-dimensional equations (|l.ip - (jl.2p with 
M/ ) = P°( a > 1/2) was studied by Li, Li and Xin in [22] and the phenomena of vacuum 
vanishing and blow-up of solutions were found there. The global existence of weak solutions 
for the initial-boundary-value problems for spherically symmetric compressible Navier-Stokes 
equations with density-dependent viscosity was proved by Guo, Jiu and Xin in [ID]. More 
recently, there are some results on global existence of weak solutions to the Cauchy problem 
(ll.ip - (ll.3p . The existence and uniqueness of global strong solutions to the compressible Navier- 
Stokes equations (|l.ip - (jl.3p were obtained by Mellet and Vasseur [31] where no vacuum is 
permitted in the initial density and for < a < ^. However, the a priori estimates obtained in 
[31j depend on the time interval thus do not give the time-asymptotic behavior of the solutions. 
The first result about the time-asymptotic behavior of the solutions to the Cauchy problem 
(|l.ip - (jl.3p is obtained by Jiu- Xin [19], where the global existence and large time-asymptotic 
behavior of the weak solutions were considered in the case that p + = p_ > and u + = u_ = 0. 
In the case that p+ = p_ > and tt+ = u_ = 0, the vanishing of the vacuum and the blow- 
up phenomena of the weak solutions were also obtained in [19] . One of the key elements in 
the analysis in [19] is an interesting entropy estimate which was observed first in [20] for the 
one-dimensional case and later established in [H [2j [3] for more general and multi-dimensional 
cases due to the structure that the viscosity coefficients vanish at the vacuum. This entropy 
estimate provides higher regularity of the density and played a crucial role in [TDJ CDS [22] for 
global existence and large time asymptotic behaviors of weak solutions. 

The stability of rarefaction waves of the ID compressible Navier-Stokes equations with 
density-dependent viscosity was studied in [18] under general initial perturbations such that 
the initial data and the solutions may contain the vacuum. However, in [18], the rarefaction 
wave itself is away from the vacuum. In this paper, we are concerned with the the case when 
the rarefaction wave is permitted to be connected to vacuum. 

For definiteness, we consider the case of a 2-rarefaction wave such that p„ = 0, p+ > 
in (jl.3p . Similar to [18], we will first construct a class of approximate solutions satisfying 
some uniform estimates and furthermore prove the global existence of weak solutions for the 
Cauchy problem of (|l.l|) - (|1.3j) . To get the uniform energy and entropy estimates in time to the 
approximate solutions, we combine the elementary energy estimates and the entropy estimates 
in an elaborate way. Note that the elementary energy estimates and the entropy estimates 
are coupled to each other due to the underlying rarefaction wave. This is quite different from 
the previous works on the global existence and the time-asymptotic behavior of the solutions 
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to Navier-Stokes equations (jl.ip with density-dependent viscosity where the elementary energy 
estimates and the entropy estimates can be derived independently. Moreover, compared with the 
case of non- vacuum rarefaction waves in [p3] , some new difficulties occur due to the degeneracies 
at the vacuum states in the 2-rarefaction wave. To overcome these difficulties, we first cut off the 
2-rarefaction wave with vacuum along the rarefaction wave curve and then derive some uniform 
estimates with respect to both the approximations and the cut-off process. More precisely, for 
any v > 0, a suitably small parameter, the cut-off 2-rarefaction wave will connect the state 
(p, u) = (v,u u ) and (p + ,u + ) where u v can be obtained explicitly and uniquely by the definition 
of the 2-rarefaction wave curve. For any fixed v > 0, one can obtain a weak solution to the 
compressible Navier-Stokes equations (jl.ip - (|1.3p with (p_,m_) replaced by (i/,uu u ) along the 
same line as in our previous paper [18J. Thus, in order to get the solution to the original problem 
(jl.ip - (|1.3p . we will derive some uniform estimates with respect to both the approximations and 
the cut-off process. To this end, the approximation parameters e and the cut-off parameter 
v should be chosen in an appropriate way. Thus, as a limit of this approximate solution, a 
global weak solution to (|l.ip - (|f -3p is shows to exist with the uniform- in-time estimates ()2.f 81) 
and (12T91) . 

Next, we study the large-time asymptotic behavior of any weak solutions to (jf .1 p - (jl.3p 
under the uniform- in-time bounds (|2.18p and (|2.19p . It is shown that time-asymptotically, the 
density function tends to the rarefaction wave connected to the vacuum in L°° norm. This time- 
asymptotic behavior of the density function implies that the vacuum in the far field is essential 
and will maintain for all the time. This is quite different from the previous results in [18] and 
|22| where all the possible vacuum states will vanish in finite time. At last, we prove that such 
a weak solution becomes regular away from the vacuum region of the rarefaction wave by using 
the Di Giorgi-Moser iteration and higher order energy estimates. 

Notations. Throughout this paper, positive generic constants are denoted by c and C, which 
are independent of e, v and T, without confusion, and C(-) stands for some generic constant(s) 
depending only on the quantity listed in the parenthesis. For function spaces, L p (£l), 1 < p < oo, 
denote the usual Lebesgue spaces on Q C R := (— oo, oo). W k,p (£l) denotes the k th order Sobolev 
space, H k (Q) := W k ' 2 (n). 



2 Preliminaries and Main Results 

In this section we first describe the rarefaction wave connected to the vacuum to the compressible 
Euler system (jl.4p . Then an approximate rarefaction wave will be constructed through the 
Burger's equation and the main results of the paper will be given at last. 



2.1 Rarefaction waves 



The Euler system (|1.4j) is a strictly hyperbolic one for p > whose characteristic fields are both 
genuinely nonlinear, that is, in the equivalent system 
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the Jacobi matrix 

u p 

y(p)/p u 

has two distinct eigenvalues 

\i(p,u) =u - vV(p), A 2 (yO,n) =u+y / p'{p) 
with corresponding right eigenvectors 



r t (p,u) = (l,(-iyy^-) t , z = 1,2, 
P 

such that 

n{ P ,u) ■ V {P)U) \i{ P ,u) = (-I)' p/(P) ^ (P) ¥= o, i = 1,2, ifp>0. 

Define the i— Riemann invariant (i = 1,2) by 

£ i (p, u) = u + (-iy+ 1 J P ^ids, 

such that 

V( Pi „)Si(yO,n) • r*(p,u) = 0, Vp > 0, u. 

There are two families of rarefaction waves to the Euler system (jl.4p - (|1.5j) . Here we consider 
the case of 2— rarefaction wave connected with vacuum, that is p_ = m_ = 0, p+ > 0. Thus 
we can define the velocity at the positive far field 144. = First we give the description of 
the 2-rarefaction wave connected with vacuum, see also in details in [25]. From the fact that 
2— Riemann invariant is constant: 

£ 2 (j0- = 0,it_) = S 2 (p+,ti+), 

we can define the velocity u_ which is the speed of the gas coming into the vacuum region. 
Then the entropy condition A 2 (p_ = 0, it_) < A 2 (/o + ,ti + ) is always satisfied. This 2— rarefaction 
wave connecting the vacuum p_ = to (p + ,u + ) is the self-similar solution (p r ,u r )(^), (£ = |) 
of JOD-drgi defined by 



u. 



p r (0 = 0, if £<a 2 (o, u _; 

£, if u- <£< A 2 (p+,u+), (2-1) 

A 2 (/9+,u+), if £>A 2 (p+,u+), 



A 2 (/f (£)X(0) = 

and 

E2(p r (0.« r (0) = S 2 (p+,n+) = E 2 (0,u_). (2.2) 
Thus we can define the momentum of 2-rarefaction wave by 

_ Jp r (£K(£), if P ^)>o, 

m (f ) = < (2.3) 

[0, if //(£) = o. 

In this paper, we consider the time-asymptotic behavior toward such rarefaction waves of solu- 
tions to the compressible Navier-Stokes equations fll.lj) with density-dependent viscosities. 
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2.2 Approximate rarefaction waves 

Consider the Riemann problem for the inviscid Burgers equation: 

wt + ww x = 0, t > 0, x G R, 

w(x, 0) = 




(2.4) 



If W- < then the Riemann problem (|2.4p admits a rarefaction wave solution «; r (x, i) = w r (j ) 
given by 



»'(f) 



10. 



f <w_ 



(2.5) 



> 



t w +- 

Consider the solution to the following Cauchy problem for Burgers equation 
wt + ww x = 0, t > 0, x € R, 
w(0, x) := wq(x) = — '—^ h 



W+ + W- W + — W- 



/ (l + y 2 )-^y. 



(2.6) 



o 



2 2 

Here q > 2 is some fixed constant, and K q is a constant such that / (1 + y 2 )~~ q dy = 1, and 

Jo 

r/ is a small positive constant to be determined later. It is easy to see that the solution to this 
problem is given by 

w(t,x) = w (x (t,x)), x = x (t,x) + w (x (t,x))t. (2.7) 
Then the following properties hold (see |30|). 

Lemma 2.1 Let w_ < w+, the Cauchy problem (|2.6p has a unique smooth solution w(t,x) 
satisfying 

i) < w(t,x) < w + , w x (t,x) > 0; 

ii) For any p (1 < p < oo), there exists a constant C pq such that 



\w(t,-) -w r (-)\\ LP{R) <C p 5 r i] p, 



Wx{t) ||lp(r)< C pq mm{5 r n l p, 6?t 1+p }, 



where 5 r 



II Wxx{t) ||lp(r)< C pq m.m{5 r r] p, r/ ^ 
w + — w_, and C p ,C pq are independent oft; 



2pq 



Hi) sup \w(t,x) — w r (j)\ — > 0, as t — > oo. 



We now turn to rarefaction waves to the Euler system f)l .4|> - fj 1 . 5 [) . Set \2(p±,u±) = w± with 
p_ = in (12, 4h . Then the unique solution (p r ,u r )(^) in (I2.1j) - (|2.3p to the Riemann problem 
(|1.4p - (|1.5p can also be expressed in terms of w r (j ) in (|2.5p . by 



A 2 (/( i ),^(|))=^( i ), 
S 2 (p r (j),n r (|)) = S 2 ( P± ,u ± ). 



(2.8) 
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Correspondingly, an approximate 2— rarefaction wave (p,u)(t,x) can be denned by 

\ 2 (p(t,x),u(t,x)) = w(l +t,x), 

(2.9) 

^2(p(t,x),u(t,x)) = T, 2 (p±,u±). 
It can be checked that (p,u)(t,x) also satisfies the Euler system 

p t + (pu) x = 

(2.10) 

(pu) t + (pu +p{p)) x = 0, 
and properties listed in the following Lemma. 

Lemma 2.2 The approximate 2— rarefaction wave (p,u)(t,x) defined in (12.9P satisfies 

7-3 _ 

(1) p x > 0, U x > 0, U x = y/jp 2 p x ; 

(2) For any p (1 < p < oo), there exists a constant C pq such that 

_i 

Up, «)(*»•) - (/^X^IIz^r.) < c P (w+ -«;_)?7 p, 

IK(V)l|i>(R) < C pg m.m{5r] p,Sp(l + t) p}, 
||« m (V)||li>(R) ^C^m^ 

where 5 = | p+ — p_ | + | u + — u- \ is the strength of the rarefaction wave; 

(3) lim sup \p(t,x) - p r (f)\ =0. 



Remark 2.1 For any 1 < p < +oo, 







(*, -)\\LP(R)dt < C, 



where C is independent ofT. Note that in the case p = 1, the constant C in the above estimates 
is not uniform in T. Moreover, the following estimate holds: 



[ \\u xx (t,-)\\ L o a(R) dt<Cr j ^ [ (l + t) 1 wdt<Cri 
Jo Jo 

where C is independent ofT. 



2 

4 9 +l 



2.3 Main Results 

Set 



p 



' p 7 _ p l _ 7p 7 l {p- p) 



(T-l)pL' 



(2.11) 
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The initial values are assumed to satisfy: 



Po > 0; too = a.e. on {x G R\po(x) = 0}; 

{pTh* e f 2 (R), po*(po,p±) e with P _ = 0; 

Po P± Po P± 



(2.12) 



where u_ := — — since p 



to_ = 0. Note that (|2.12p implies that po G (7b (R) which is the 
space of bounded and continuous functions. 

Equivalently, the assumptions (|2.12p can be rewritten as 



Po > 0; too = a.e. on {x G R|po(^) = 0}; 
(Po _l )xG L 2 (K), p *(po,Po)eL 1 (R); 



(2.13) 



p (— -n ) 2 G^(R), po^-^gi^R) 
Po Po 



,m 



where (po, uq) '■= (p, n)(0, x) is the initial values of the approximate 2— rarefaction wave (p, u)(t, x) 
constructed in (I2.9j) . 

Before stating the main results, we give the definition of weak solutions to (|l.ljl - (jl.3p with 
the far fields (p±,m±). Let T > be given. For any far fields (p±,u±) satisfying p± > and 
any smooth functions (p,u)(t,x) connecting with them, we define 



Definition 2.1 A pair (p,u) is said to be a weak solution to (ll.ll )-( fl~3l) with the far fields 
(p±,u±) provided that there exists a smooth functions (p, u)(t, x) with the same far fields (p±,u±) 
and p > 0, such that 

(1) p > a.e., and 

p€ L oo (0,T;L oo (R)))nC([0,r] x R), pV(p,p) G L°°(0, T; L X (R)), 
{p a ~K G F°°(0,T;L 2 (R)), x/p(«-«) G L°°(0, T; L 2 (R)); 

f^J For any £2 > *i > an d £OT 2/ C £ Co(R x [£1, ^2])? mass equation (jl.ip /10/ds in t/ie 
following sense: 



I p(dx\\\ = I [ (p( t + y/py/pu( x )dxdt; 
JR Jti JR 



(2.14) 



("3; For any V G C^R x (0,T)), ii holds that 



Jo Jn {^^ U ^ + + 



dxdi + (p a u x ,ip x 



(2.15) 



where the diffusion term makes sense if written as 



{p a u x , ip x ) = - J J p a 2 ^fpuijj xx dxdt - - J" - J J (p a 2) x ^/putfj x dxdt. (2.16) 



The first main result in this paper reads as 
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Theorem 2.1 Let a and 7 satisfy that 

1<7<2, and 1 < a < , (2.17) 

and suppose that (|2.13j) holds. Then the Cauchy problem (|l.lj )- (|1.3|) admits a global weak solution 
(p(x , t) , u(x , t)) in the sense of Definition \2.1\ with the smooth function (p,u) replaced by the 
approximate rarefaction wave (p,u). Furthermore, this weak solution (p(x , t) , u(x , t)) satisfies 

p>0, max p<C, p € C(R x [0, XI), (2.18) 

(i,t)6Rx[0,T] 



sup 

te[o,T] 



[IVp(«-«)| 2 + (p a -^ + p*(p,p3]dx 

7+a 1 7+q i (2-19) 

+ / / {[{p 1 ^ - p 2± ^~ A ) x } 2 + u x p^(p,p)+pu x (u-u) 2 + A(x,t) 2 }dxdt <C, 
Jo Jr 

where C is an absolute constant depending on the initial data but independent ofT and A(x,t) 6 
L 2 (R x (0,T)) satisfies 



p^ Aipdxdxt = — / / p a ~2 y/p(u — u)(p x dxdt 



Jr Jo Jn, 

T 



2a 



[ [ {p a ~^)xVp{u-u)ipdxdt, V93 G Co°(R- x (0, T)). (2.20) 
Jo Jr 



2a - 1 7 y R 

Remark 2.2 should be noted that there is no requirement on the sizes of the strength of the 
rarefaction wave and the perturbations. The class of initial perturbations given by (|2.13|) is quite 
large compared with those for the case of constant viscosities, I26\j , I29f . \30). 



Remark 2.3 The important case of the shallow water model, i.e., a = 1,7 = 2, is included in 
our theorem. 



Remark 2.4 For any far fields (p±,u±) satisfying p± > and any smooth functions (p,u)(t,x) 
connecting with them, one can also obtain the existence of weak solutions in the sense of Defi- 
nition \2~l\ in a similar way (see 123)/ for the weak solutions in the case of p± > and u± = 0). 
However, in order to get the uniform in time estimates in ()2.18p -( i2.19p . it seems that the far 
fields (p±,u±) should be specified and in Theorem \2.1\ the smooth function (p,u)(t,x) is replaced 
by the approximate rarefaction wave (p, u) . 

The next result concerns on the asymptotic behaviors of the weak solution, which can be 
stated as 



Theorem 2.2 Let a andj satisfy (|2. 17j) and suppose that (|2.12p holds. Suppose that (p,u)(x,t) 
is a global weak solution of the Cauchy problem (ll.ip -( fl~3l) in the sense of Definition \2. 1\ satisfying 
(I2TT8]) and fl2Jgp . Then it holds that 



Consequently, 



lim sup \p(t, x) - p(t, x)\ = 0. (2.21) 

t->+oo xgR 



lim sup|p(x,t)-p r (f)| = 0. (2.22) 
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Remark 2.5 A direct consequence of Lemma \2.2\ (2), (|2.19p and (|2.22p is the following time- 
asymptotic behavior of the density function: 

lim |K-»*)-/(t)ILlp(R)=0, V2<p<+oo. (2.23) 



Remark 2.6 Theorem \2.2\ implies that for any weak solution (p, u) to the Cauchy problem 
(jl,ip - (jl.3p with the far fields given by the vacuum state and (p+,u+), if (p,u) satisfies the 
bounds (|2.18p and (|2.19p . i/ien i/ie density function converges to the 2-raref action wave to the 
corresponding Euler equations connecting the vacuum state and (p + ,u + ) in sup-norm as t tends 
to infinity. Consequently, the initial vacuum at far field will remain for all the time, which 
is contrast to the case of non-vacuum rarefaction waves studied in fT8\j where all the possible 
vacuum states will vanish. 



Finally, we can obtain the following higher regularity to the velocity function u(t, x) to a 
global weak solution (p,u)(t,x) of the Cauchy problem (|1.1|) - ()1.3[) in the sense of Definition l2.il 
satisfying (I2.18P and (I2.19P in the region away from the vacuum region of 2-rarefaction wave 



Theorem 2.3 (Regularity of the solution away from the vacuum) Let (p,u) be a weak solution 
to the Cauchy problem (ll.ip -( fL3|) satisfying (|2.18p and (I2.19p . For any fixed a > 0, there exist 
a straight line x = \%t with A2 = M(p, ^)|(p,u)=(o-,u (T ) defined in (|5.ip and a large time T a , such 
that if (t,x) E fl a := {(£, x)\t > T a ,x > A|i}, then the density has the lower bound 

P(t, x ) > |- 

Furthermore, for any (£*,a;*) G Cl a and for any r, s > such that Q* s := B r (x*)x (t*, C il a 

with B r {x*) being the ball with the radius r and the center x*, there exists a constant oq E (0, 1), 
such that 

u E C£°- (Q* s ), u E L% C (U,U + s,HUB r (x*))), 
ut E Lf oc (Q* ), u E L 2 Xoc {U,U + s,H? oc (B r (x*))). 



3 Existence of a weak solution 

We first study the following approximate system: 

Pt + {pu) x = 0, x E R, t > 0, ^ 

{pu) t + (pu 2 +p{p)) x = {Pe(p)ux)x, 

where p £ (p) = p a + ep e with e > and 9 = ^. This kind of the approximation ep e was first used 
in [T7] and 9 = \ is crucial in getting the lower bound of the approximation density function to 
enure the existence of approximate solutions. 

To overcome the difficulty caused by the vacuum in the rarefaction wave, we first cut off the 
rarefaction wave along the wave curve. More precisely, for any v > suitably small and to be 
determined, let (y, u{v)) be the state such that 



T, 2 (u,u(u)) = S 2 (/3+,n+), 



11 



where 2-Riemann invariant Ti2(p,u) = u — ^jP 2 • Then (v,u{v)) is connected to (p+,u+) by 
a non- vacuum 2-rarefaction wave given by (p r v ,u r u ){j). Then it holds that 

\{p r M)(j) - (p r v ,m r v ){j)\ < Cu. 

One can compute that u(y) = ^j^ 2 ^ - + E2(p+, «+)• So the corresponding smooth approximate 
2— rarefaction wave (p v ,u v ) described in Section 2.2 can be constructed by setting 

\ 2 {v,u{v)) = w-, X 2 (p+,u + ) = w+. 

Consequently, (p u ,u u ) will converge to (p,u) point-wisely as v tends to zero. In fact, v = v{e) 
will be chosen suitably such that v(e) — > as e tends to zero. 

The initial values (po,mo) can be regularized in a similar way as in [18] such that 

{(u, i/u(u)), as x — > — oo, 
(3.2) 
(p+, m + ), as x — > +oo, 

and 

1 

PoeA x ) > min{i/, ^e 3 — 1 }, Vx G R, (3.3) 

for suitably small e, v > 0. 
Furthermore, poe,^ satisfies 

Poe,i>*(po e ,i/, Poi/) -> Po^PcPo) in ^(R), (Po £ TJ /2 )^ ^ {Pq~ 1/2 )x in L 2 (R). 

2 

Since in the following paper a,-j satisfy (|2.17p and we choose v = £3, it holds that 



£'[{\np Qe>v ) x Ydx = { T Y / p^[(p, e J) x Ydx<C. 

R Ot — 2 JR 



While moe.j/ satisfies 



POeA Uqu) -> f>o{ u Q ) m L (R), 

POe^ PO 



and 

P0e.^( 

POe.i/ PO 



PoeM i*ok) ->-Po( «o) m L (R). 



For any fixed T > and £^v > 0, we will first construct smooth approximate solutions 
(p £tU ,u £tU )(x,t) to (|3.ip with initial values (p, pu)(0,x) = (po e ,u,moe,u)(x) defined in [0,T]. To 
do this, a key step is to get the lower bound of the density. Then the global existence of 
weak solutions to (|l.ip - (ll.3p can be proved by compactness arguments. We intend to deduce 
the uniform energy and entropy estimates with respect to e, v such that one can pass to the 
limit £,v — > 0. Due to the closeness to the vacuum of the rarefaction wave, we will have 
to combine suitably the elementary energy estimates with the entropy estimates to get the 
following estimates which are crucial to prove our main results. 
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Lemma 3.1 Suppose that the conditions in Theorem \2.1\ are satisfied and that (p £ ,u,Ue,u) is a, 
smooth solution to ([3.ip satisfying p e ^ u > 0. Then for any T > and e, v > with e ln(l+T) < C 
and e, v sufficiently small, the following estimate holds 



\pe,u(u £ ,u ~ u u ) 2 + (pe,v 2 ) + e 2 [(lnp etU ) x ] 2 + p E)U y(p Ejl/ ,p u ) \{x,t)dx 
te[0,T] JR 1 L J J 

+ / / \{Uu)x P{Pe,u) ~ P{pu) ~ P{pu){Pe,u ~ Pu) + Pe,u(u u )x(u £ ,u ~ U u f (3.4) 
JO JR I L 



{u e , v - U u ) x + (p, 



-y + CK — 1 7 + a — 1 



P 



V )x 



(x, t)dxdt < C. 



where C > is a universal constant independent of e, v and T. 

In the following, the subscripts e,v in the approximate solution (p £tU ,u e)U )(t,x) and the 
subscripts v in the approximate solution (p u ,u u )(t, x) will be omitted for simplicity. 
Proof: Step 1. Energy Equality 
It follows from (|3. 1 |) q that 

pu t + puu x + p(p) x = {p e {p)u x )x- (3.5) 
Subtracting f|3.5j) from the second equation of (|2.10p gives 

p(u - u) t + pu(u - U) x + (p(p) -p(p)) x + (P~ P)ut + {PU ~ PU)U X = (p £ (p)(u - u) x ) x + (p e (p)u x ) x . 

(3.6) 

Multiplying f|3.6j) by u — u yields 



p(u — n) 2 i v pu{u — u) 



+ 



+ (U - u)(p(p) - p(p))x ~ Pe (P) (u ~U){U- U) 



+Pe(p) (u - U) a 



= p £ (p) U x (u - U) + p £ (p)u xx (u - u) - (p - p)u t + (pU - pu)u x (u - u) . 
J x L J 

(3.7) 

Note that ^(p,p) defined in (|2.1ip satisfies 



P(p)x 



pV(p,p) + puV(p,p) +{u-u) x {p{p)-p{p))+u x p 1 - p 1 -7P 7 l {p- p) 



(p- p)(u-u). 



(3.8) 



P 

It follows from ((321) and flM} that 
' p(u — u) 2 , , .1 „ / 

2 + P*(P, P)\ f + X) + p e (p) (U - 6) a 



P 1 -f ~ IP 1 (P ~ P) 
P(P)a 



p £ (p) U x (u - u) + p £ (p)u xx (u - U) - (p - p)Ut + (pU - pu)u x -\ — (p-p)(u-u), 

P J 

(3.9) 



where 



pu(u — u)^ 

Hl(t,x) = h pu^(p,p) + (U- U)ip(p) -p(p)) - Pe(p)(u ~ U)(u - U) x . 



Since 



(p - p)u t + (pu - pu)u x + (p - p) = p(u - u)u x , 
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we obtain 

~ p(u — u) 2 s i . , .r, 

2 + P*(P, P)\ t + Hi*(t, x) + ii e {p) \{u - u) x 

+p(u-u) 2 u x = p £ (p) u x (u-u) + p £ (p)u xx (u-u). 



Step 2. Entropy Equality 
Rewrite (13.61) as 



+ u x 



(3.10) 



p(u-u) t + pu(u-u) x + (p{p)-p(p)) x + (p- p)u t + (pu- pu)u x = (p a l + ep e l )pu x . (3.11) 

J x 

Note that 



where ip £ ,0 (p) with 9 = | is defined by 



-p(^' e (p)U-pu(^ e (p)U, 



(3.12) 



cpf(p) 



< p^ 1 p e 

T^i +£ T- 



l 



In p + e 



P 



, if a/l,Q>0, 



if a = 1. 



Thus (|3.1ip becomes 



p(u - u) t + pu(u - u) x + (p(p) -p(p)) x + (p- p)ut + (pit - pu)u x = -p(<pf 6 (p)) xt - pu{(p^ £ {p)) xx . 

(3-13) 

Multiplying (I3.13P by (ife' {p)) x shows that 



P{^\p)) 



+ 



pu(ipf' 6 \p))l 



+ 



p{u-u)^f\p)) x + pu{u-u){^/{p)\ 



-(u-u) [p(^ e (p)) xt + pu(^ e (p)) xx \ + (^ e (p)) x (p(p) - p(p)) x 
+(ip^ e (p)) x (p- p)u t + (pu- pu)u x = 0. 

Combining (j3TT3l) with (ETUI) yields 

[\p [(« - u) + (^(p))*f } t + {±pu [(« - n) + 

+ (.<Pe' e (.P))x{p(.P) ~ P{P))x + {U-U) (p- p)u t + {pu - pu)u x 
+ (<PT e (p))x (p ~ p)u t + (pu - pu)u 



(3.14) 



} +(u-u){p{p) -p(p)) a 

J X 



0. 



(3.15) 



Step 3. A Priori Estimates 

It follows from (gTS]) and fl3~T5]) that 



(u - u) + (<p? e (f)) x +pu^(p,p) 



[\p [(« - ii) + (^(p))*] 2 + p*(p, p)} f + { \pu 
+(u - u){p{p) - p{p))\ + u x p{p) - p(p) - p'{p){p -p)+ p(u - ufu x 

J x L J 

+( l fT (p))x (p - P)u t + {pu - pu)u x + p(p) x - p(p) x = 0. 



(3.16) 
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Now we deal with the last term on the left hand side of (|3.16p . Note that 



(p - p)u t + (pu - pu)u x + p(p) x - p(p) x = p(u - u)u x + p{p) 3 



pp(p% 



and 
Thus 



J -2 



(<Pe' e (p))x (p ~ PH + (PU - pu)u x + p(p) x - p(p) x 



-+4 

a 



(u - u)u x + (p a 2 p x + ep y 2 p x ) p(p) x 



0-2 



PP(p)a 



(3.17) 
(3.18) 

(3.19) 



Direct computations show 
PP(p)x 



P a ~ 2 Px 



P{p)x 



P 



, J 2 - 2 

(a + 7- l) 2 

ZPj , a+7— 1 , 7 — a — 1 

(p 2 )xP 2 



1 2 



+ 



O'y a + 7 — 1 a+7— 1 a+7 — 1 . 

(p 2 )x(p 2 -P 2 



(a + 7- 1) 



a(a + 7 — 1) 
a(a + 7 — 1) 



(P a " P a 



a+7 — 1 -7 — a — 1 

(P^~~ )xp 2 



x (a + 7 — l) 2 
(P a -P°). 



a+7— 1 a+7— 1 a+7— 1 

(P 2 )«e(p 2 -P 2 ) 



(3.20) 



and 



P 6 ~ 2 P* 



P(p)x 

"47 



pp(p)a 



+ 



9+7-1 fl+7-1 
2 — o 2 

+ 7-I) 2 

27 9+7-1 _ 7-e-i e 

(p 2 )xP 2 {p -F) 



7 , 9+7-1 , . 9+7-1 9+7-1 , 

(P 2 )x{p 2 -P 2 , 



+ 



0(0 + 7-1) 
2 7 



+ 7-I) 2 

87 



+ 7-1)' 



, 9+7-1 

r(P 2 



9+7-1 9+7-1 . 
.(pa -p 2 ^ 



(p 



9+7-1 7-9-1 
2 )zP 2 



(/-A 



0(0 + 7-1) 
Substituting (l3~T9l) - (l3T2!l into (l3~T6|) gives 

{^P [(« " u) + &e' e (p))x] 2 + p^(p, p)} t + H 2x (t, x) + u x [p(p) - p{p) - p'(p)(p - p) 

47 



(3.21) 



+p(u - ufu x + ( — + e^— ) (u-u)u 
\ a r 



+e 
+e 



47 


(0 + 7- 


I) 2 


87 




I) 2 






0(0 + 7" 


-1) 



9+7-1 9+7-1 



(P 

9+7-1 
(P 2 



P 



(a + 7- l) 2 L 

87 , a+- 

(P 2 



a+7— 1 a+7— 1 . 

(p 2 - p 2 _ 



)xx{p 

V— 1 7-9-1 



(a + 7- l) 2 

7-1 9+7-1 , 27 
2 — p 2 



-(P 



a + 7— 1 a+7— 1 



a(a + 7 — 1) . 



a+7— 1 7— a— 1 

(P 2 )xP 2 



) 

(P a -P a ) 



(3.22) 
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where 



H2(t,x) = -pu 

+ 



(u-u) + {<p? e (p)) x \ + pu*(p, p) + {u- u)(p(p) - p(p)) 

O'y a + 7— 1 a + 7 — 1 a + 7 — 1 

(P 2 )x(p 2 -P 2 ) 



(a + 7- l) 2 

y^'V Q+7 — 1 7 — a — 1 

(p^-^P^^-p-) 



+e 



a{a + 7 — 1) 
87 



(3.23) 



(e + 7-1) 2 
27 



, fl+7-l . . 9+7"! fl+7-1 . 

(p 2 MP 2 -P 2 . 



y— 1 7-9-1 



Multiplying (|3.22p by a and then adding up to (|3.10p and noticing that p e {p) = {p a )x+£{p e )x 

L J x 

in the right hand side of ()3.10p . one can get 



a 



{u-u) + {<pf{p)) x 



2 p(u - u) 2 

+ = + 



(a + l)p*(p,p)| + aH 2 (t,x) + Hi(t,x) 



+(a + p(p) - p(p) - p'(p)(p - p) + (a + l)p(« - u) 2 ^ + (p a + ep 9 ) (it - u) a 



+ 



4aj 



(q + 7- If 



Q+7— 1 Q+7 — 1 , 

(P 2 -P 2 



+ e- 



4cry 



+ 7" I) 2 



9+7-1 9+7-1 . 
(P 2 -P 2 



a . 



+ 



p a u xx (u - u) +eyp u xx (u -u) + (1 - -)(p ) x (u - u)u x 

807 , a+7-1 , a+7-1 n+7-1 807 
(P 2 )xx(p 2 -P 2 , 



(q + 7- l) 2 

27 

(a + 7- 1) 



Q + 7 — 1 7 — a — 1 

[P 2 )xp 2 



(p*-p*)-e- 



+ 7-I) 2 
2cry 

? + 7-l) 



9+7-1 9+7-1 9+7-1 . 
(P 2 )«e(p 2 -P 2 . 



, 9+7-1 . 7-9-1 
(P 2 )*P 2 



(p* 



Integrating (j3.24p over [0, t] x R with respect to t, x gives 

2 p(u-uf 



IJ: 



a 

1 «• 2 



(u-u) + (^(p)) a 



+ 



(a + l)ptf(p,p)}(t,x)da 




+ 

/o ^/r 

+(p a + ep; 
4ct7 



(a + l)it x p(p) - p(p) - p'{p){p - p) + (a + l)p(it - u) 2 ^ 



u — u) 



+ 



4cry 



(a + 7 - 1)' 



q+7— 1 Q + 7— 1 . 

(P 2 -P 2 . 



+£ 



(0 + 7 — l) 2 



9+7-1 9+7-1 , 
(P 2 -P 2 . 



a 

2 P0 



(uo-u ) + (tp?> e (p )), 



jdxdr 

2 po(-up - up) 2 



(a + l)p ^(po,Po)}^ 



+ 1, 



(3.24) 



(3.25) 



where 



p a u xx (u -u)+e p u xx (u - u) + (1 - -)(p ) x (it - u)u 



+ 



°8q7 



(a + 7- l) 2 
27 

(a + 7- 1) 

6 

:=£'- 

8=1 



(/> 



Q+7— 1 



Q + 7 — 1 Q + 7 — 1 

.(pa -p 2 ) +e - 



807 



Q+7— 1 7— Q— 1 

{P~^)xP~ 2 ~ 



(P a -P a )~e- 



7 + 7-1) 



, 9+7-1 9+7-1 9+7-1 . 

(p 2 Mp 2 -p 2 . 



. 9+7-1 , 7-9-1 

(P 2 )*P 2 



(p e -p e )}dxdr 



(3.26) 
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We now estimate the right hand side of (|3.26j) terms by terms. First, 



h= / p a u xx (u — u)dxdr 
lo Jr. 



= / \/p(u — u)p a 2u xx dxd,T 

Jo.Jr^ (3.27) 

ny/p(u - u)p a ~^u xx [l\ {0 < p < 2p+} + l\ {p > 2p+} ]dxdT 
i 

■ = hi + h2, 

where and in the sequel l|n denotes the characteristic function of a set C (0, t) x R. 
Rewrite I\ 2 as 



'12 




^p(u - u)u xx [(p a *-p a i) + p a 2}l\ {p > 2p+} dxdT 
o Jn (3.28) 



T 1 4- T 2 

M2 T J 12- 



Using Lemma 12.21 (and its Remark 2.1) and noting that a > 1, one has 

~ u )\\L 2 (~R)\\ u xx\\L 2 (~R)dT 

JO 



<C sup \\y/p{u - n)|| L 2 (R) / \\u xx \\ L 2(u\dT (3.29) 
te[o,T] Jo 



< C SUp ||Vp(«-M)lli»(R), 

te[o,T] 



and 



< C SUp ||VP( n -^)llL2(R) SUp \\(p a '^ - P a ~^)l|{p>2p+}I|L2(R) • / \\Uxxh°°(R)dT 

te[o,T] *e[o,T] jo 

<C??35+t SU p ||^(-u-n)|| L2(R) sup ||(p Q -5 -p a -5)l| {p > 2p+} || i2(R) 

te[o,T] te[o,T] 

<Ci]^ \ sup \\Vp{u-u)\\ 2 L2(R) + sup -^-5)11^ } ||2 

L te[o,T] *e[o,r] 

(3.30) 

Note that if a and 7 satisfy 

7 + 1 

1 < a < -L^-, (3.31) 

then 2(q — i) < 7, and then 



1 *x-h\2 ~ 1 



Km ^ _ ijm (7 - 1)( ^, p> < c (3 32) 

P^\P, P) p^+co pi — pi — jpl L (p — p) 



Thus if 1 < a < -^-i— , then 



sup ||(p a -2-^-2)l| {p >2p + }lli 2 ( R )<C sup ||p*(p,p)|| Ll(R)! (3.33) 
te[o,T] te[o,T] 



for some uniform constant C > 0. 
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Substituting ([338]) . ([339]) . (ICTjl and (f3T33l) into (I33TD yields 



h<Cr)^ sup ||v^(«-«)ll£2 (R) + sup ||p*(p,p)|| £ i(R) 
L te[o,T] te[o,T] 



+ C„. 



(3.34) 



Next, Io can be rewritten as 



h =e 




10 JR 

:= hi + -^22 



First, since 9 = i, it follows that 



p 9 u xx {u - u) + (1 - -^){p e ) x {u - 

'0 JR 

-p e u xx (u -u)-(l- -^)p 9 (u - u) x u 3 




dxdr 
dxdr 



(3.35) 



hi =£ / 2au xxy /p(u - u)dxdr 

JO JR 

< Ce sup - «)IIl 2 (R) / ||«xx||i,a(R)dT 

iefo,Tl Jo 



On the other hand, 



'22 



te[o,T] 

< Ce sup \\y/p{u - u)\\ L 2 {n) . 

te[o,T] 



<7 T / /[(u-ny 2 dxdr + Ce f / 
4 Jo Jr jo jr 



(3.36) 



/3 u x dxdr, 



(3.37) 



while 



/ / p 8 u 2 x dxdT=ef I p e [l\{o<p<2 P+ } + l\{ P >2 P+ }]u 2 x dxdT 
Jo Jr Jo t jr 

= £ Io J 4^ ei|{ °- p - 2p+}+ [^~^) + ^] 1 l^>2p+}}^^ 

< Ce\u{l + T) + Ce sup ||(/ - pP)l\ {p > 2p+} \\ LH n) [ ll^lll-rR,)^ 
te[o,T] Jo 



<Ce\u{l + T) + Ce sup ||p*(^p)|| L i (R) , 

te[o,T] 

where in the last inequality we have used the fact that 

p^+oo p^S>(p,p) 



(3.38) 



since 



~ < 1 < 7. 



Substituting the estimations (|3.36|) - f[3.38j) into (|3.35p . one can get 

I 2 <CeHl + T) + - I I P e [{u-u) x ] 2 dxdT + Ce sup \\ P V(p, p)|| Ll(R) 

4 JO JR *G[0,T1 



+Ce sup 11^(1/ — ix)||x,2 (B .). 
te[o,T] 



(3.39) 



It follows from the fact that 



g+7— 1 g + 7 — 1 27 

n 2 — O 2 a+7— 1 

lim ^- £ ^ = 1, 



p^o p&(p,p) 
that for any e > 0, there exists 5 e > 0, such that if < p < 5 e , then 

Q+7— 1 a + 7 — 1 2 7 

I l/o 2 ~P 2 T+ 7 ' 1 _ n < e 
P^(p,P~) 

Fix e = i, then there exists <5i > 0, such that if < p < Si , then 

z 2 2 

a + 7— 1 a + 7— 1 2 7 

|/0 2 — p 2 | a+7-1 1 

' P*0M) 11 < 2' 

thus for any p > 0, 

^P*(p,P~) < Ip 2 ^ -p'^I^T < -p^(p,p), if 0<p<<5i. 
Similarly, it follows from the fact that 

Q+-7— 1 a 4- 7 — 1 ^7 

p 2 — p 2 U+7-1 

lim —, — rr = 1, 

P^O pV(p,p) 

that there exists 81 > 0, such that if < p < 61 , then 

2 2 

a + 7 — 1 a + 7 — 1 2 7 

|p 2 — p 2 J a+7-1 1 

1 P*0M) 11 < 2' 

thus one can choose ^ < <5i such that for any p > 0, 

1 _ g + 7— 1 a + 7 — 1 27 3 _ - 

-p^{p,p) <\p 2 -P 2 h+T" 1 < -p*(p,p), if f<p<<fi. 

The term ^3 can be estimated as follows. Since 

, a+7-1 . .a + 7 — 1 a+7-3 a + 7 — l,_a 

P 2 )«s = ( o P 2 = — o F (P 2 ^ ^ 

2 2^/7 

a + 7 — 1 _a _ a(a + 7 — 1) a+1-7 _ 9 

= P 2 ^xx + P 2 <, 

2^7 4 7 

and 

a+l-7>2-7>0, 

one can rewrite ^3 as 

f l f / 4a,/7 _a_ 8a 2 a+ i_ 7 \ a + 7 _i a+7-1 

13 = / / — 1 — 7P 2u xx + — : T Kp 2 Up 2 -pa )cbdr 

JoJr v " + 7-1 a + 7-1 / 

<C/ f \{u xx ,u x )\\{p^^ - p^^^dxdr 

v<P<5i} 



t 



= C / \(u xx ,U x )\\(p 2 -p 2 )|(l| {0 <p<<h} + 1 \{8 1 <p<2p + , 
JO JR 27 

+ 1 l{<5i<p<2p+, 5i</5<p+} + 1 l{p>2p + })^T 
:= ^31 + -^32 + -^33 + -^34- 
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Direct computations lead to 



hi < C I I 11(^,^)11 27 1 l )l\{o<p<8 l }\\ 27 dxdr 



g-f-"7~ 1 g+'y — 1 



<C sup 2 -p 2 )l|{o<p<5 1 }ll 2^ / [1(^,^)11 2 7 dr 

|[3T43l "+7-i 

< C sup ||p*(p,p)|| 2 - 



(3.48) 



te[o,T] 
a + 1 



-— sup \\pV(p,p)\\ L l(VL)+C a . 

8 te[o,T] 
Similarly, due to f|3.46j) . one has 

/ 32 <^ sup \\ P *(p,p)\\LHn) + C a . (3.49) 
8 te[o,T] 

On the other hand, 

/33<C/ / ||(^,^)|| L2(R) ||( p ^-_p^-)l| {5 

JO JR 2 2 

Q!-|-'y — 1 (X-\-'~f — 1 / rt 

<Csup \\(p 2 -p 2 )l| {(5l <„<2p +) Ji^p+llli^R) / ll(^,^)||L2(R)dr (3.50) 
£6[0,T] 2 2 jo 

<C sup ||p*(p,p)||| 1(R) , 
te[o,T] v ' 

where one has used the fact that 

n+7-i Q+7-1 , n i 
\P 2 _ P 2 J 1|{<Ji <p<2p+, 5i <p<p+} 
2 1 < (J 

P^(P,P) 

Moreover, I34 can be estimated as 

hi < C [ [ ||(w xx ,^)|| LOO(R) ||(p £; ^ -p £ ^)l {p > 2p+} || L i (R) (ixdT 



JR 



rt 

OL-\-'~f — 1 Ct-|-'Y — 1 / rt 

< C sup IK/O 5 -p-^~ )1 {p >2p + }||li(r) / 11(^^,^)11^°° (R.)^ (3-51) 
te[o,T] Jo 

<Cr]^ sup ||p*(p,p)||£i(R), 
te[o,T] 

due to the facts that 

a + 7 — 1 a-f-'Y — 1 

hm l^-^l 1 ^} < ^ 
p->+oo p^(p,p) 

since a < implies 

< 7, i.e., a<7 + l. 

Now we turn to the term 1$. First, 

a+7-1 7-a-l ] Q! + 7 — 1 , 7-1 , 

(P 2 )xP 2 = r- \P 2 U x ) x 

\x 2^/7 

2^/7 47 v ; 
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Thus 




2l 



aH-7 — 1 7 — a — 1 

(p-^—) x p-^— ( p Q - p a )dxdT 



< 



R (a + 7- 1) L 

C / / 1(^,^)11^° - P a )|(l|{0<p<2p+} + M{p>2p+})dxdT 
JO JR. 



'■= hi + -^52- 



One has 



hi < C f \\{u xx ,u 2 x )\\ L 2 m \\{p a - p a )l\{o<p<2 P+ }\\mn.)dxdT 
t 

<C sup \\{p a -P q )1|{o< p <2 P+ }IIl2(r) / IK^,^)|| i2(R) cZT 
te[o,Tl J 



< C sup ||p#(p,p)||| 1(R) , 
te[o,T] w 

due to the fact that 

( P a - p a ) 2 l\ {0 < p <2 P+} _ (7-l)(p a -p") 2 l| { 0<p<2p +} 

^ Ap 2a - 2 + P 2a - 2 )(p-p)\ , 

- ° Tp _ pf i l{0<P<2p+} 

<C, 

where one has used 1 < 7 < 2 and a > 1. 
On the other hand, 

h2<C [ \\(u xx ,ul)\\ L oo {R) \\(p a -P a )|l|{p>2p + }||L1(R)^T 
JO t 

<C SUp ||(P Q - P a )l|{p>2p+}||L1(R) / ll(^x,^)|| L oo (R) dT 

te[o,T] jo 
^Cr/^r S up ||p#(p,p)|| L i (R) , 

*6[0,T] 

due to the fact that 

lim \p a - p a \n {P >2 P+} = lim (7-i)iP"-P a iii { p> 2 p +} 

p->+oo P^(p, p) p-++oo p7 — pi — 7pT _1 (p — p) 

< C, if a < 7. 
Finally, we estimate the terms I4 and Iq. As for ^3, one has 

, 9+7-1 , .0 + 7 — 1 9+7-3 _ ^ + 7 — 1,9 

(P 2 )xx={ P 2 Pzjx = - .- (P*U X ) X 

2 2^7 
6* + 7-l_9_ 6>(6> + 7 - 1) e+1-7 2 

= - P 2 ^xx + "j P 2 Uj, 

2^/7 47 
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therefore, 

h<Ce[ [ \{u xx ,p^^ui)\\{p^^ - p^^)\dxdt 
Jo Jr. t 

< Ce^V / / \{u X x,ul)\\(p~3 J1 - p~~ a^^dxdt 



-i 



= Ceu 2 / / \(u xx ,u x )\\(p 2 - p 2 )|[l| {0 <p<| } + l|{f<p<2 P+ } + l\{p>2 P+ }]dxdt 

JO J R. 

:= Iai + ^42 + Iaz- 

(3.52) 

Recall the following useful fact that for any given C > 0, there exists a constant /3 G [0, 1] such 
that 

P^(P,P)l|{0<p<C} = P J ~ P 1 -IP^ip- P) ^{O^KC} 

= [pp + {l- p)pf-\p-pfl\ {0 ^ c} (3-53) 
> max{p+, C} 7 ~ 2 (p - p) 2 , 

provided that 1 < 7 < 2. 
One can compute that 

U\<Cev-^ I I Uuxxiul^Ltm^lip 2 ^ -p^ _L )l{o< P <^}llL2(R)^r 

<Cei/V S up \\(p^r- -p^r- )1 { o< p <^}||l2(r) / life, ul)\\ L 2 iR) dr (3.54) 
tefo.Tl jo 



^C^K - 1 ) sup ||p*(p,p)||| 1(RV 
te[o,T] y ' 

where one has used the facts that 

fl+7-l 9+7-1 -9 , , 9+7-1 9+7-1 ,9 , 

(p 2 _p 2 ri {0 < p <| } (p 2 - p 2 ri {0 < p< | } 

p*(p,p) - (p + )-^(p-p)2 ' 1 I ' <1 ' 

and the function 

fl+7-l 8+7-1 ,9 

(p 2 -p 2 Y 



1 



(P-P) 2 

is monotone decreasing in p £ [0, |], that is, 

9+7-1 _ 9+7-l -o 9+7-1 9+7-1 „ 

2 If 2 ) W|l < Um (p^~P^) 2 = , +7 -3 < ^-1 
(p - p) 2 p^O (p - pf 

Moreover, it holds that 



I 4 2<Ceis 2 1 f f ||fe,x4)|| L 2 (R )||(p '2 -p 2 )l{|<p< 2p+ }||£a(R)da?rfT 

9-1 9+7-1 8+7-1 /"* _ _„ 

<Cev 2 SU p || (p 2 - p 2 )1 { -< p <2 p+ }||l2(R) / IKttaas.uJUia^jdr 
«e[o,T] Jo 

KCeuW-V sup ||p-p|| L 2 (R) 
te[o,T] 

(£55) 3 „ „ I 
< CW^D sup ||p*(p,p)||| 1( 



(3.55) 
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where in the third inequality one has used the fact that there exists a constant /3 S [0, 1], 



(p 



9 + 7-1 9+7-1 , 



L {|<P<2p+} 



(p - py 



= (e + 1 -if[p P + {i-p)p] 
<{e + i-i)%) 6 - 1 . 



M-y-3. 



L {|<P<2p+} 



And then 



h 3 <CEV^ I / || ( p - 

'o Jn 



y-1 9+7-1 



P 2 ) 1 {p>2p + }\\L^R)\\{Uxx,U x )\\ L oo {n) dxdT 



9-1 9+7-1 9+7-1 / _ _„ 

<Cei/2 sup ||(/3 2 -p 2 )1 {p>2p+ }||li(r) / ||(^,%)IU°°(R)rfr 
te[o,T] Jo 

<Cev-^ sup ||p*(p,p)|Ui(R), 
te[o,T] 



(3.56) 



since 



9+7-1 9+7-1 , , 

hm (P 2 'j/ >>^>=Q. 



In summary, by combining (|3.54p . (|3,55p and (|3.56p . one can arrive at 



sup ||p*(p,p)||li( R ) + 1 

te[o,T] 



(3.57) 



Finally, Iq can be estimated similarly as for I± and the details will be omitted for brevity. 

For definiteness, we take ev~^~^ = ea, i.e., v = £3 since Q = \. Consequently, choosing e 
such that eln(l + T) < 1 and e suitably small and combining all the above estimates shows that 
for a and 7 satisfying (I2.17p . 



+ £■ 



[ (in p) J 2 +p*(p, p)}(x,t)dx 



sup / \p(u - u) 2 + \(- r 

te[o,T] Jr l lv " - 2 

/ {«x P(p)-P(p)-P'(p)(p- P) + p(u - u) 2 u x + {p a + ep e ) (u-u)^ 



JR 



(3.58) 



a+7— 1 a+'y- 

{p~^~ - p~ 



) z | +e\(p 2 -p 2 ) X J 



todt < C. 



Thus Lemma 13. II is proved. □ 
The following lemma is the key point to get the existence of the approximate solution 

(Pe,v,U e ,v)(t,x) with V = £3. 



Lemma 3.2 There exist an absolutely constant C and a positive constant C(e, u, T) depending 
on e, v and T such that 



0<C(e,u,T) < p £yU <C. 



(3.59) 



Proof: From the Gagliardo-Nirenberg inequality: 



L°°(R) < CH/xlliafRjII/llip^,)) 
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where < /3 < 1, 1 < p < oo to be determined, and /?, p satisfy 

P _ 1-/3 
2 p ' 

we have 

SUp -p Q "5|| ioo(R) 

te[o,T] 

<c sup \\(p a -^ - P a -h,C (R) sup iip a -^-^ni; ( % 



< Cjup [||(p a -^),||J 2(R) + IK^-^xll^^J jup ||p a "* -p* 



i£[0,T] 



< C SUp ||/9 a 2 — p" 



te[o,T] 
due to the fact that 



i - 1 II 

Hlp(R)' 



2 

Since 



»)*IIl3(r) = !!(«- 2u x || L 2 (R) < C\\u x \\ L 2 {R) < C. 



(3.60) 



_ 1 _ 1 27 

\o a 2 — n a 2 I 2o-l 

lim £ -f-^J = 1, (3.61) 

p^O p#(p,p) 



there exists a positive constant 5\ , such that if < p < 5\ , then 

2 2 
1 1 2 7 

I- (3-62) 

p^{p,p) 2 

thus, for any p > 0, 

ip*(p,p)< Ip^-l-^-^l^ <^*(p,p), if < p < . (3.63) 
2 2 2 

Similarly, there exists a positive constant 5\ , such that for any p > 0, 

2 

Jp*(p,p) < -ff~k\2§h < lp*(p,p), ii 0<p<d\. (3.64) 

2 2 2 

Set 

Then for such p = and 7 6 (1,2], one has 

p*(p,p) 1 K5\<P<2p + , 5\<P< P+ } < C 1 kn<P<2p + , 6\<p< P+} < C, (3.65) 

and 

lim . r 1 = lim ^ ^ — ' =7-1. 3.66 

p->+oo pW(p, p) p^+00 p~t — p~? — p~f L [p — p) 

Thus it follows from ([3.66P that 

|p«-3-p«-3|f <Cp*(p,p), ifp>2p + . (3.67) 
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Collecting (13T63D . (IOTD . (1051) and (13^71) gives that 



R 



/ |P° 2 -P° TCMftK/K^} + % 1 1 <p<2p+, 0<p<*\} 

■/a 2 4 4 (3.68) 

+ 1 l{5 1 1 <P<2p+, 5i<p<p+} + 1 l{p>2 / 9+}) da; 



2 2 



< C f p^(p,p)dx < C. 
Jn 



Substituting (|3.68p into (|3.60p yields the uniform upper bound for p £>v (t,x). 

Next we derive a lower bound for p £jU {t,x). Since lim^o p^f(p,p) = ff > v 1 ', then p^>(p,p) 
is bounded away from on [0, ^p\. Thus one can deduce from the bound on pty(p,p) in 
L°°(0, T; L 1 (R)) that there exists a constant C\ = C\{v,T) > 0, such that for all t G [0, T], 

1 _ If 

meas{x G R|p(x,t) < -p(M)} < — — -- / pV{p, p)(x,t)dx < C x . 

2 mt pe[o,|p] /^W 9 ' ^{xeR| P (x,t)<|p(x,t)} 

Therefore, for every xq G R, there exists M = M(u, T) > large enough, such that 

/ Pe,u(x,t)dx > / p £>u (x,t)dx 

J\x-x \<m J{\x-x \<M}n{xen\p s ,„(x,t)>±p(x,t)} 

> - inf t)meas< {\x — xq\ < M} n {x G R|p E £) > -p(x, t)} > 
2 (x,t) I 2 J 

= ~meas|{|a; - z | < M} n {x G R|p e ,„(x, t) < -p(x, i)} c } 

> ^(2M-d) >0, 

for all t G [0,T]. 

Due to the continuity of p £jV , there exists x\ G [xq — M, xo + M] such that 
PeA x i,t) = TH Pe,u(x^) dx > 7T7( 2M ~ 

IM J\ x - X0 \<m ^ M 

Thus, 

\\n.p e>t/ (x ,t)\ = \]np EjV (xi,t)+ / (\np £tV ) x (x,t)dx\ 

J X\ 

< I ]np S)V (xi,t)\ + ||(lnp £ ^) x (-,i)|| L 2 (R) |xi - x |^ 

< C(e,u,M,T) +C £ M^. 

Consequently, we can get that there exists a positive constant C(e,u,T) such that 

Pe,u(x ,t) > C(E,V,T), 

for any xq G R and t G [0, T]. □ 

With the lower and upper bounds on p £)U , we can get the existence of the approximate 

solution (p £u ,u £U )(t,x) by a similar argument as in [32]. In order to pass the limit e — > with 

2 ' ' 

v = £3, we need the following higher estimates on the momentum. 
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Lemma 3.3 There exists a positive constant C(T) independent ofe,is, such that 

sup / p e ,u\ u e,v ~ u w \ 3 (t,x)dx +/ / p E: „[(u £jU - u„) x ] 2 \u £>w - u\dxdt < C(T). 
te[o,T]JR J o JR 

The proof of Lemma 13.31 can be done along the same line as in our previous paper [18] and we 
omit the details for brevity. 

Now with these uniform in e, v estimates at hand, we can pass the limit process e — > 
with ^ = £3, obtain the existence of the weak solution (p,u)(t,x), and get the uniform in time 
estimates in Theorem 12.11 



4 Asymptotic behavior of weak solutions 

In this section, we will study the asymptotic behavior of a given weak solution (p,u)(t,x) to 
the Cauchy problem (jTTjl - lfOl) in the sense of Definition O satisfying (pTTH]) and (pTHJj) . We 
assume that the solution is smooth enough. The rigorous proof can be obtained by using the 
usual regularization procedure. 

Proof of Theorem 12. 2i For any s > 1, by the uniform upper bound of p, it holds that 

\p s -p s \ 2 <C\p-p\ 2 . 

Hence it follows from (j2~T9j) and (j333|) that 

f \p s - p s \ 2 dx <C [ (p- pfdx <C. (4.1) 

JR JR 

Similarly, 

f \p s - p s \ 2X dx < C [ \p- p\ 2X dx <C [ (p- pfdx <C, (4.2) 
Jn Jn Jn 

for any A > 1. 

Set b = 2±2=1. Then one gets from (I2J9D and (13331) that 

f I W ~ P%? + MP - P?}dxdr < C. 
Jo Jn 

For s > b + 1 , it holds that 

(p s -p s )\t,x)= f [{p s -p s )\dx 

J — oo 

=2 r & - - ?) x dx 

J — oo 

= 2, f (p s - - p 6 )^" 1 + {p b ) x {p s ~ b - p s ~ b )]dx 

J — oo 

<C||p s -p s ||L 2 (R)ll(p fe - A'IIl 2( r)+C / p?u x {ft-f?)(jf- b -p 

Jn 

<\\P S ~ p S \\LHR)\\{p b ~ p b )x\\LHK)+C / U x (p-p) 

Jn 



)dx 
dx. 



Consequently, 



f sup(/) s - p s ) A dt 
Jo ieR 



< c sup i|p s -/5 s ni 2(R) 



f\\{p b 

fe[o,T] v ' J o 

+C sup ||p-p||i 2(R) / / 
te[o,T] Jo JR 



P )z|Il 2 (R) 



u x (/0 — p) 2 dxdr 



< C. 



Moreover, applying (j4.2|) leads to 



(p s - p s )\p s - p s ) 2l dxdt 



< 



sup (,0 s -p s 



\2I 



R 



< sup f (p s - p s ) 2l dx [ sup(p s - p s fdt <C, VZ > 1. 
te[o,T] jr Jo xeR 



Set 



Then 



/(*) = / (p s - p s ) 4+2 ^x. 

JR 

/(*) G L 1 (0,oo)nL oo (0,oo) 



due to and flOp . 

Furthermore, direct calculations show that 



j t f{t) = (4 + 2/> - p s f +2l (p s - l p t - f-^dx 

= -(4 + 2l)s I (p s -pn 3+2 V-Hpu) x -p s -Hpu) x }dx 
R /• 

= (4 + 20(3 + 2l)s / (p s - p s ) 2+2l (p s - p%{p s - l pu - p^p^dx 

+ (4 + 2Z) S ( S -l) / {p s -p s f +2 \p s - 2 PxP u-p s - 2 p x pu)dx 
. JR 



(4 + 2/)(3 + 2Z)s / {p s -p s ) 2+2l (p s -p s ) x p s {u-u)dx 



R 



+(4 + 2Z)(3 + 2Z) S / (p s -p s ) 2 + 2l (p s -p%u(p s -p s )dx 

+(4 + 2Z)a(s-l) / (p s - p s f+ 2l p s - l p x {u - u)dx 

+(4 + 2Z) S ( S -l) / {p s -p s f +2l u(p s - l Px -p s - 1 p x )dx 
Jn 

■= Ji(t) + J 2 (t) + J 3 (t) + Ji(t). 
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Now we 



claim that Ji(t) E L 2 (0, +00), (i = 1,2,3,4). In fact, 

^ (4 + 20(3 + 20^ f ( )2+2 y (u __ } r s - 6(6 _^ + ( s - b _ - s - b){ ^] dx 

b Jr L 

f4 + 2/)f3 + 2ns 2 /■ , . ..0,0, _. 



2 Jr 

= (4 + 20( b 3 + 2 ° g jtf - P s ) 2+21 Mu ~ u)p^ b -H P b - P%dx 
(4 + 20(3 + 20s 2 r , . 



IV 

/ (P s - P s ? +21 Vp(u - u)p s ^(p s - b - p s ~ b )(p h ) x dx 
JR 

\\(n b - n b )J\ r o,„, 



Thus, 



JR 

<C\\^p(u-u)\\ L 2 {R) \\(p b -p b ) x \\ L2m 

+C\\^P(U - U)\\ L 2 {R) \\(P S - p S ) 2+2l ( P S ~ b - P^Xp'UmR) 

< C\\y/p(u - u)\\ L 2 {R) \\(p b - p b )x\\ L 2(R) + C\\y/p(u - u)\\ L 2 {R) \\u x (p - p)|| L 2 (R) . 

(4.6) 

> b ) X II 2 T2lT>\dt 



f\Ji(t)\ 2 dt<C sup \\^-p{u-u)\\l 2{R) f \\{p b - p b ) x \\h(R) 
Jo te[o,T] v 'Jo 

ft 



+C sup \\Vp( u ~ u )\\h(R) [ |K(P-P)H|2(R) 
te[o,T] Jo 



dt 



UP* ~ P%\\h(K)dt 



< C sup ||^(u-u)|| 2 2(R) / 

te[o,T] Jo 

+C sup \\<fp(u- n)||| 2(R) / / u x (p- pfdxdt 
te[o,T] Jo Jr 

< C. 



The fact that Jj(i) G L 2 (0,+oo), (i = 2,3,4) can be shown similarly. 
Thus 

|/(t)GL 2 (0,+oo). 
Combining the fact that f(t) G L x (0, +00) n L°°(0, +00), one has 

/(i) ->• 0, t ->• +00. 
be determined later. One has 



Let m > 1 be any real number to 

\(p s -p s r\ = \ f [(p s - p s r%dx\ 
j —00 

(p s -p s r~\p s -p%dx\ 



= \m 
= \ms 



dx\ 



[ x (p s - p s r~ l \ — -—i(p a ~^)xp s ~ a+ ^ - p 3 - 1 ^ 

J — 00 2 

II(/9 Q ~^)x||l2(R) + \\P~* L Ux\\v>(EL) 



<c\\( P s -p s r- i \\ LHR) 

<C\\(p s -p s r- X \\ L2{R y 



Choosing 2(m — 1) = 4 + 21 yields 

-P S T\ < Cf^(t) -> 0, as t^+oo. 



sup I Go* -p s 

x£R 



(4.7) 



(4.8) 



(4.9) 



(4.10) 
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Therefore, 



lim sup \p s — p s \ = 0. 



Now we prove that lim sup \p — p\ = 0. 
Since 

lim ~. ~ — ^j— - = 1, uniformly in p, 

p->o+ |p s — p s \ 

then for any a > 0, there exists 5 a > 0, such that if < p < 5 a , then 

i \P~P\ S n . 

i — t — 1 < a. 

\p s — p s \ 

Thus, fix a = i, then if < p < 5 := 5i , one has for any p > 0, 



\P-P\ s <\\p s -P s \. (4.11) 



Now 



\P - PY =\P~ Pl'(l{0<p<5} + 1 {p>S,0<p<§} + 1 {p>5,p>|}) 
3 

< " ^fe/K*} + Cs\p s - P s |l { p>5,o< P <f } + C«5lP s - P s \ Sl {p>5, P >^y 

Therefore, 

sup |p - p\ s < Cg sup \p s - p s \ + Cs sup \p s - p s \ s -> 0, 
xGR i£R i6R 



as t — )■ +oo, which implies that 

The proof of Theorem 12.21 is finished. 



lim sup \ p — p\ = 0. 



5 Regularity of the solution away from the vacuum 

In this section, we will prove Theorem I2.3[ that is, we will show that away from the vacuum 
region of the 2-rarefaction wave (p r , u r )(^), any weak solution (p, u)(t, x) to the Cauchy problem 
(|l,ip - (jl.3p satisfying ()2. 18[) and (|2.19p becomes regular as stated in Theorem 12.31 

Due to the definition of the 2-rarefaction wave p r (t,x) in (12.8H . for any fixed a > 0, there 
exist a unique u a such that (p, u) = (a,u a ) lies on the 2-rarefaction wave curve. In fact, 

2-v/t 7-1 

u a = — —v 2 +Z 2 (p+,u+), (5.1) 
7-1 

since 2-Riemann invariant ^(p, u) = u — ^y/? 2 " 5- is constant along the 2-rarefaction wave 
curve. Due to the expanding property of the 2-rarefaction wave, we have the lower bound of 
the density function p r (t,x) of 2-rarefaction wave on the right of the straight line x = \%t with 
= A 2 (cr,n (T ) = A 2 (p,u)| (([VU ) = ( . U(r ), that is, 



p r (t,x)>a, if x>\ a 2 t. 



(5.2) 
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Thus it follows from the asymptotic behavior (|2.22p of p(t, x) that for ^ > 0, there exists a 
large time T a such that if t > T a , then 

||p(V)-/(V)|U~(R) < \- (5.3) 

Therefore, in the domain 

n a := {(t, x)\t > T a , x > Aft}, (5.4) 

it holds that 

p(t,x) > -, if (t,x) £ n a . 

So in the domain Q a , any vacuum states vanish and thus the higher regularity of the weak 
solution (p,u)(t,x) can be expected as stated in Theorem 12.31 In the following, we give the 
proof of Theorem 12.31 First we establish the local uniform boundedness of the velocity u(x, t) 
by the De Giorgi-Moser iteration method. To this end, we rewrite the momentum equation as 

u t + uu x + 7P 7_2 /0 X = p a ~ 1 u xx + ap a ~ 2 p x u x , (5.5) 

For any (t*,x*) G fi<j> and for any r, s > such that Q* s := B r {x*) x + s] C fi CT , and 

for any test function C( x >t) 6 ^^{Qrs) satisfying < £ < 1, one can get from the uniform 
estimates in (pHS]) and (pTT9|) that 

sup / (u-u) 2 +p 2 x + (p-p) 2 (x,t)dx+ f I (u-u)l + (p-p) 2 x dxdt < C. (5.6) 

t,<t<t,+s J B r (x,) L Jt, JB r (x,) L 



It follows from the construction of the rarefaction wave (p,u)(x,t) that 

sup I {u 2 + p 2 x + p 2 ){x,t)dx + [ [ {u 2 x + p 2 x )dxdt < C '. (5.7) 

t,<t<t,+s J B r (x,) Jt, JB r (x,) 

Multiplying the equation (|5.5[) by C 2 (u — k) + for any fc£R and integrating the resulted equation 
over B r (a;* ) x (t* , t] for i £ (i* , + s] , one arives 



1 



C 2 (^~ A;)^(x,t)dx+ / / p"" 1 ^- k) + ] 2 x dxdt 

B r {x,) Jt* JBr(x t ) 

If C 2 (u-k) 2 + {x,U)dx+ ( ( {cQtiu - k)% + l(( x [(u - k). 

1 JB r (x,) Jt, JB r (x,) K 6 



<B r {x,) Jt, JB r (x„) 

- !P^- 2 PxC 2 (u - k)+ + p a ^C 2 (u - k) 2 + + p a - 2 p x u x Q 2 (u - k) + }dxdt 
<\( . C\u-k) 2 + (x,Qdx + l [ [ p a - l [t{u-k) + ] 2 x dxdt 



B r (x,) ° Jt, JB r {x,) 



(5.8) 



+c f I {(\Ct\ + \C x \ 2 )(u-k)l + C\( x \[(u-k) + } 3 + \ Px \c 2 (u-k) + 

Jt, J B r {x,) k 

+\ Px \ 2 ( 2 (u- k) 2 + }dxdt, 
where in the last inequality one has used the fact 
\p a ~ 2 p x u x ( 2 (u - fe) + | < C\p x \\(u x \((u - k )+ 

= C\p x \\((u) x - ( x u\((u - k) + 

< C\p x \\(( U ) x \((u - k)+ + C\p x \(\( x \(u - k) 2 + 

= C\p x \\[C{u - k) + ] x \t{u - k)+ + C\p x \C\Cx\(u - k)\ 

< \p a - l [Q{u - k) + ] 2 x + C\p x \ 2 ( 2 (u - k) 2 + + C\Cx\ 2 {u - k)\. 
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Thus from (|5.8p . it holds that 

/ C 2 (u-k) 2 + (x,t)dx+ f [ 

J B r (x,) Jt, JB 



B r (x,) 



[C(u — k) + ] 2 dxdt 



< C 



B r (x,) 



C 2 ( U -k) 2 + (x,Qdx + C / \(\Ct\ + \Cx\ 2 )(u-k) 



(5.9) 



t, J B r (x,) 



+C\Q[{u - k) + f + \p x \( 2 (u - k)+ + \p x \ 2 Q 2 {u - k)%}dxdt. 
Now the last three terms in the last integral of (|5.9p can be estimated by 




'U J B r {x, 




H, JB r (x,) 

< c 

< C 

< C 

< C 

< C 



C\Cx\[{u - k) + } 6 dxdt < C / \\Q x {u - k) + \\ L 2 (Br(xt)) \\((u - k)i\\ L 2 {Br{Xt)) dt 



\C x \ 2 {u-kf + dxdt + C 



It, J B r ( X , 




't, JB r ( X , 



It, J B r ( X , 




t, J Bf( X , 



t* JB r (x,) 



C 2 (u-k)idxdt 



\( x \ 2 (u-k) 2 + dxdt + C \\((u-k) 



u 



\L°°(B r {x*))\ 



( u - k )+\\h(B r (x*)) dt 



\( x \ 2 (u-k) z + dxdt + C / IIC^-^+lll-^^))^ 



\Q {u-k) + dxdt + C J \\((u- k) + \\ L 2 {Br{Xt)) \\[C(u- k) + ] x \\ L 2 (Br(Xii)) dt 
\( x \ 2 (u-kY + dxdt + ^ [ \\l((u - k) + } x \\ 2 L2{Br(xt)) dt 



+C \\((u-k) + \\l 2{BriXf)) dt, 



(5.10) 



t, J B r (x,) 



\p x \C 2 (u - k) + dxdt 



\ \ IAz|C 2 (^ - k) + dxdt 

J JQ* >3 n[u>k] 

< C||C(u- k) + \\ L6{Q , s) \\p x \\ e 



(5.11) 



L5 (Q* s n[u>k]) 



< 



C\\C(u - k) + \\ V2{Q , a) \\ Px \\ L 2 {Q , a) \Q% n[u> k]\s 



and 



t, J B r (x,) 



< -||C(n - k) + \\ 2 V2{Q , e) + C\Q% n [u > # , 
\Px\ 2 C 2 (u - k)\dxdt 



< 



C r SUp||p a .|| i a( Bp ( a! ,)) / ||C(«- k )+h^(B r (x,)) dt 



(5.12) 



<C \\((u- k)+ \\ L 2(B r (x,)) II [C(« ~ k )+\x \\tf{B r { x ,))dt 



< 



\\[C(u-k) + } x \\i HBr(Xf)) dt + C / K(u-Q + \\UBr(x.)) 



sdt. 



Note that in (IBTTTT) . the space V 2 (Q* jS ) is defined by V 2 (Q* >S ) = {/ € L 2 (Q* rs ) \ \\f\\v 2 (Q* s ) < 
+oo}, with the norm 

11/11 V 2 (Q* J =ess sup ||/(-,<)||L3(B t .( a! .)) + ll/*IU a (Q;.)- 
te[t,,t,+s] 
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Substituting (i57L0l) - ([5J2l) into (\5JB implies 



[ ( 2 (u- k) 2 + (x,t)dx + [ I [((u - k) + ]ldxdt 
<[ ( 2 (u- k)\{x,U)dx + C [ [ ( 2 (u-k) 2 + dxdt (5.13) 

J B r (x f ) JU J B r {x*) 



+i||C(«-fc)+llv 3W ..) + C / / ^ (\Ct\ + \Cxi z )(u-k)ldxdt + c\Q:. s n[u>kp. 

Applying Gronwall's inequality to f|5. 13j) gives 



sup [ ( 2 (u - k) 2 + (x,t)dx + f [ 

U<t<U+s J B r (x*) JU JB r (x*) 

rt*+s 



[C(u — k) + ] 2 dxdt 



<[ ( 2 {u-k) 2 + (x,U)dx + C [* [ (\Ct\ + \Cx\ 2 )(u-k) 2 + dxdt + C\Q%n[u> fc]|§. 

JB r (x*) JU JBr(x,) 

(5.14) 

With the estimate (I5.14h at hand, one can show that u is bounded above locally by the clas- 
sical De Giorgi-Moser iteration method and choosing suitably k. Similarly, one can obtain the 
estimates to (u — £;)_ as in (|5.14|) . Thus we can get the lower bound for u locally. Furthermore, 
one can get the local Hoder estimates of u by the classical parabolic theory, that is, there exists 
a positive constant ao £ (0, 1), such that 

u€Cg'^(0*,). 

In the following, we will further show that the weak solution u(x,t) is in fact a strong solution 
locally as stated in Theorem 12.31 Rewrite the momentum equation as 

pu t + puu x + "/p~'~ 1 p x = p a u xx + ap a ~ x p x u x . (5.15) 

Multiplying (|5.15|) by Q 2 ut and integrating the result over B r (x*) x (i*,i] with t G (t*,i* + s), 
one can get 

p( 2 u 2 dxdt 

U J B r (x,) 

= 11 ( - pu( 2 u t u x - jp' y ~ 1 p x C 2 ut + p a ( 2 u t u xx + ap a ~ 1 p x u x ( 2 u t )dxdt 

Ju JB r {x,) v ^ ' (5.16) 

< \ f I P C 2 u 2 t dxdt + C f [ (pu 2 ( 2 u 2 x + p^-Zplt 2 

1 Ju JB r (x,) JU JB r {x,) v 

+p 2a - l C 2 u 2 xx + p 2a -\ 2 p 2 x u 2 x )dxdt. 

By the uniform upper bound and lower bound of the density p(x, t) in the domain f2 CT and the 
local boundedness of the velocity u(x,t), it follows from (|5.16p that 

f [ C 2 u 2 dxdt <C f [ U 2 u 2 x + Cpl + C 2 u 2 xx + C 2 p 2 x u 2 x ) dxdt 

JU J B r (x t ) Ju J B r (x*) v 7 

<C I j \c 2 ( u -u) 2 x + C 2 (p-p)l + C 2 (u 2 x + p 2 x )]dxdt (5.17) 

JU JBr(x,\ L J 



't* J B r (x t ) 

7 



-elf 
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Note that 



it Ib ^ p2 ^ dxdt - j t \\P*\\h(B r (x*))\\C-Ux\\ 2 L ™ {Br (x*)) dt 



SUP \\Px\{i*(B r (x*)) / \\(Ux\\L*(B r M)\\(CUx)x\\L*(B r (x*))dt 
t&(t,,t,+s] 

I' 1 



[ [ C 2 u 2 xx dxdt + C p f [ (C 2 u 2 x + ( 2 u 2 x )dxdt 

Jt, J B r (x,) Jt* J B r (x,) 

n( 2 u xx dxdt + Cp, 
-l r (x,) 



(5.18) 

where f3 is a small positive constant to be determined and Cp is the positive constant depending 
on p. 

Thus it follwos from (J5TTD , (foTTH]) with /3 = 1, and (pUT)) that 

( 2 u 2 t dxdt<C + C [ [ Q 2 u 2 xx dxdt. (5.19) 

It, J B r {x,) Jt, J B r (x,) 

Multiplying (|5.5j) by C 2u xx and integrating over B r (x*) x (t*,t] with t G (i*,i* + one can get 




p Q 1 ( 2 u xx dxdt 

't, JB r (x,) 

C 2 u t u xx + ( 2 uu x it xx + 7/? 7 ~ 2 ''Px( 2 u X x - ap a ~ 2 p x u x Q 2 u xx )dxdt. 




t, JB r {x,) 



(5.20) 



Note that 



( 2 u t u xx = (( 2 u t u x ) x - ( 2 u x u xt - 2(( x u t % 

\ 
2 



C 2 u 2 ( 5 - 21 ) 
(C 2 u t u x ) x - (— -^) t + (Ctul - 2(( x u t u x , 



thus it holds that 

"t r /-2 n ,2 



[ [ ( 2 u t u xx dxdt = - I ^-^(x,t)dx+ [ [ (CCtu 2 x -2CCxU t u x )dxdt. (5.22) 

Jt, Jb t (x,) JB r (x,) 2 J tm JB r ix,) 



>t, JB r (x,) JB r (x,) z Jt, JB r {x,) 

Substituting (j5T2"2"j) into (^2"UD gives 



\ [ ( 2 u 2 x (x,t)dx + / / ( 2 ul x dxdt 

z JB r (x,) Jt, JB r (x,) 



r(x,) Jt, J B r (x,) 

alCCtul - 2(( x it t u x + Q 2 uu x u xx + 2p x ( 2 u xx - —u x C, 2 u xx )dxdt 
Jr(x,) V P 

-22 , /-2„,2 



dxdt. 



<P (CX + ( 2 u 2 xx )dxdt + C P / (Ct + Q)< + C(P Z X + <) + c>> 

Jt, Jb t {x,) Jt, Jb t {x,) l 

(5.23) 

Combining the estimates (|5.18p . ()5. 19j) and ()5.23|) and choosing both f3 suitably small, we can 
get 

If ( 2 u 2 x (x,t)dx+ [ [ Q 2 {u 2 + u 2 xx )dxdt<C. (5.24) 

z JB r {x,) Jt, JB r {x,) 

This completes the proof of Theorem 12.31 
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